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Abstract 

We study the supersymmetric circular Wilson loops in A/" = 4 Yang-Mills theory. 
Their vacuum expectation values are computed in the parameter region that admits 
smooth bubbling geometry duals. The results are a prediction for the supergravity action 
evaluated on the bubbling geometries for Wilson loops. 
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1 Introduction and summary 

In a holographic correspondence, a quantum gravitational system is exactly equivalent 
to a non-gravitational theory on a lower dimensional space. In known examples [H [2] the 
non-gravitational systems are gauge theories. Gauge invariant observables are mapped 
to non-normalizable deformations of the gravitational background. 

Following progress in understanding the gravity duals of local operators in a gauge 
theory, the recent years have seen further extension of the dictionary to the realm of 
non-local gauge invariant operators. The picture that has emerged is strikingly universal 
for all local and non-local operators. Consider operators of the form TrR(...), where the 
trace is evaluated in a representation R of SU{N). We will often use the symbol R 
to denote the associated Young diagram as well. For different representations R, the 
operators are best described by different objects on the gravity side. The fundamental 
representation R = O corresponds to a fundamental string. A high-rank symmetric 
representation corresponds to a D-brane, while an anti-symmetric representation is given 
by another type of D-brane. Finally, an operator with a large rectangular Young diagram 
is dual to a smooth bubbling geometry with a flux supported on a new cycle. A general 
Young diagram corresponds to a combination of these objects. This pattern, summarized 
in Figure [H has been demonstrated for local operators [SI HI El E] and Wilson loops 
0[9l[IOl[Ill[l2l[T3l[Tl[T5l[T6l[I7]inA/' = 4 Yang-Mills, as well as for Wilson loops in 
Chern-Simons theory [21 [181 UHl [SO]- Higher dimensional non-local operators [211 l22l [231 
[Tni[211[2m in gauge theories have descriptions in terms of branes and bubbling geometries. 

The identification of the operators and their dual bubbling geometries has been made 
on the basis of symmetries and charges. An interesting test of the identification, and 
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D3-brane on 




Bubbling geometry 



Figure 1: Gravity duals of Wilson loops in A/" = 4 Yang-Mills. A string corresponds to the 
fundamental representation, a D3-brane to a symmetric representation, a D5-brane to an 
anti-symmetric representation, and a bubbling geometry to the representation specified 
by a rectangular Young diagram. (D3-brane D5-brane ...) in the figure should be 
replaced by (D3-brane wrapping C AdS^, D3-brane wrapping C S^) for local 
operators in A/" = 4 Yang-Mills, and by (D-brane, anti-D-brane) for Wilson loops in 
Chern-Simons theory. 



AdS/CFT itself, would be to compare the expectation value of an operator with the 
on-shell supergravity action evaluated on the bubbling geometry. For most operators, 
however, the test is trivial, the expectation value vanishes for local operators, and does 
not depend on the coupling for the straight Wilson line. Even if the expectation value does 
not vanish, the computation in the gauge theory is usually done only in weak coupling. 
It is then hard to make comparison with the strong coupling computation performed on 
the gravity side. 

In this paper we make a non-trivial prediction for the bubbling geometries dual to 
particular Wilson loops in A/" = 4 Yang-Mills, by computing the vevs of the operators 
on the gauge theory side in the strong coupling regime. The operators are the circular 
super symmetric Wilson loops defined as 

WR = TiRPexp ^{A + e%'ds). (1.1) 

Here A = A^dx'^ is the gauge field, (z = 1, 6) are the real scalars, s is a parameter 
for a circle in Euclidean such that ||(ia;/(is|| = 1, and (6^) is a constant unit vector 
in R^. A smooth bubbling geometry is expected to be dual to a circular Wilson loop 
whose Young diagram has long edges. More precisely, if we parametrize the diagram as 
in Figure [H nj for I = 1, g + 1 and kj for I = 1 , . . . , all have to be of order for 
the geometry to be smooth. The circular Wilson loop is very special in that such strong 
coupling computation from gauge theory is possible, due to the conjecture [26] that the 
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Gaussian matrix model captures the circular Wilson loops to all orders in and the 
't Hooft coupling A = Qym^- 




Figure 2: The Young diagram R, shown rotated and inverted, is specified by the lengths 
rij and kj of the edges. Equivalently, nj and kj denote the lengths of the black and white 
regions in the Maya diagram, n^+i is defined by nj = N. 



Relegating the computations to later sections, let us simply summarize the main result: 
in the limit — > oo with rii/N and kj/N kept finite, and for the large finite 't Hooft 
coupling A, the vev of the circular Wilson loop in the Yang- Mills with gauge group SU (N) 
is given by 

{Wr) = explj^'^^MKi-lRl/Nf + OiNHogX)^. (1.2) 

Here \R\ is the n-ality, i.e., the number of boxes in R. We have also defined Kj = 
kj + ki^i + ... + kg ioT 1 < I < g and Kg^i = 0. An important check of the result is that 
it is invariant under complex conjugation R ^ R of the representation, or in terms of 
the parameters, under 

ni ng+2-i, kj kg+i^j. (1.3) 

We expect that the Wilson loop vev is related to the on-shell supergravity action S'sugra 
as 

{Wr) ~ e-^="^■■^ (1.4) 

Note that {Wr) = e^^^^'^'' in our regime. The quadratic dependence on of the exponent 
is what is expected from the supergravity action S'sugra = (1/2^^^) / + ■■■ with ~ 
g^a'^, to be evaluated on the dual geometry. 

This result is obtained by studying the Gaussian matrix model that captures the 
correlation functions of circular Wilson loops. In fact we turn the Gaussian matrix 
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model with operator insertions into multi-matrix models whose partition functions are 
the Wilson loop vevs. We perform the saddle point analysis and propose the eigenvalue 
distributions. Based on this proposal, we are able to calculate the Wilson loop vevs with 
the result (Oil . 

The circular Wilson loops in Euclidean can also be thought of as operators in 
the A/" = 4 Yang- Mills defined on S"^ [26]. The loop is now the equator of 5"*, and 
the 50(2) X 50(3) subgroup of the 5*0(5) isometry group is preserved. The 50(2) is 
part of the preserved subgroup 5f/(l, 1) of the conformal group 50(1,5). The operator 
(11. ip also preserves an 50(5) subgroup of the 50(6) R-symmetry group. Thus the 
bubbling geometry duals of the Wilson loops must solve the BPS equations of the type 
IIB supergravity with the EAdS2 x 5^ x 5^ ansatz. Note that we need to work with 
Euclidean signatures since the matrix model captures the circular loop in Euclidean R^ 
only. The resulting geometries should be the Euclidean continuation of the solutions 
recently found in [T7] . 

Given the result of the present paper, the challenge is to use the dual geometries to 
reproduce the prediction (11. 2p by computing the on-shell supergravity action [27]. One 
would need to deal with such subtle issues as volume regularization, counter terms, and 
ambiguity of the action due to self-duality of the 5-form fiuxE] 

Though the analysis of the matrix models for A/" = 4 Yang-Mills is in principle sufficient, 
much confidence in the results and the proposed eigenvalue distributions comes from the 
study of surprisingly analogous matrix models that describe the Wilson loops in Chern- 
Simons theory [281 EHl [30]. For these matrix models, the resolvents and the spectral 
curves can be computed exactly in the limit oo, with rij/N and kj/N kept finite, 

without the assumption that the 't Hooft coupling is large. The eigenvalue distributions 
turn out to be qualitatively the same for the A/" = 4 Yang-Mills and Chern-Simons theory. 

The paper is organized as follows. In section [21 we review the correspondence between 
the circular Wilson loops in A/" = 4 Yang-Mills and the observables in the Gaussian 
matrix model. In section [3] we analyze the Wilson loop with a rectangular Young diagram. 
Section HI deals with the Wilson loops in general representations that admit smooth dual 
bubbling geometries. 

2 Gaussian matrix model for circular Wilson loops 

It is believed [261 ISI] that the correlation functions of circular loops in A/" = 4 Yang- 
Mills are captured by the Gaussian matrix model. The precise correspondence states in 

^ The parallel problem of matching the Wilson loop vev and the closed string partition function has 
been solved for the Chern-Simons/conifold duality to all orders in genus expansion |18| . The matching 
of the Wilson loop vevs and the dual on-shell D-brane actions has also been demonstrated [TUl HH [13] 
for A/" = 4 Yang-Mills. 
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particular that 

Tr^Pexp ^{A + ^ ~Z j ^^^^P (^X'^^^^^) '^^^^^ ' ^'^'^^ 

The left-hand side is the normalized expectation value in the Yang-Mills with gauge group 
U{N). The right-hand side is normalized by using the partition function Z which is the 
integral without the insertion of Tr^e^^. dM is the standard hermitian matrix measure. 
In the absence of operator insertions, the eigenvalues are distributed according to the 
Wigner semi-circle law in the large limit. The eigenvalue density p{m), normalized so 
that J p{m)dm = 1, is given by 

2 

p{m) = — -VA — m^. (2-6) 

ttA 

We are also interested in the SU{N) Yang- Mills theory, for which the conjecture in 
[26] states that 

Tr^jPexp ^{A + 9'X'ds)'^ =^ j dM exp (^-'^TtM^^ Tr^e^^'. (2.7) 

Note that on the right-hand side we need the traceless part 

M' ■=M- —TtM ■ In (2.8) 
of the matrix M so that e*^' G SU{N). The relation 

Tr^e^^' = (dete^O"""^'^Tr^e'^ (2-9) 

will be useful later. 

Matrix integrals have been studied intensively in the past, partly motivated by non- 
critical string theory [321 E3] and relation to supersymmetric gauge theories [Ml 155] . 
In these contexts one does not encounter operators of the form Tr^^e^^ with the Young 
diagram R that is large both in the row and column directions, and it appears that there 
is no study of them in the literature|§ 



3 Rectangular Young diagram 

In this section, we formulate 2-matrix models whose partition functions are the vev of 
the Wilson loop with a rectangular Young diagram shown in Figure [31 We will have two 
such matrix models. The first captures the geometric transition of D5-branes wrapping 
C in AdSr, x S^. The second is similarly interpreted in terms of D3-branes wrapping 
S"^ C AdS^. The techniques we introduce here are a generalization of the methods used 
in [13]. The pictures that emerge are similar to those of [121 136] . 

^ I thank E. Martinec for a discussion on this point. 
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n 



k 

Figure 3: A rectangular Young diagram with n rows and k columns. 



3.1 Eigenvalues as D5-branes 

We begin with the Yang-Mills with gauge group U{N). Let U he a unitary matrix and 
consider a sum over all the Young diagrams for which the summand is non-vanishing: 

^TRe^^TiRTU = det(l + e^' ®U). (3.10) 

R 

Here is the transpose of R. This identity, expressed in terms of Schur polynomials as 
YsR{x)siiT{y) = Yl{^ + XiVj) ioT X = (xi) and y = (yj), (3.11) 

R i,j 

is well-known [37]. The relation can be inverted: 

Tr^e^ = J dU det{l + e^^ ® f/-^)Tr^Tf/. (3.12) 

Here dU is the Haar measure on the unitary group normalized so that J dU = 1. We now 
specialize to the rectangular Young diagram R with n rows and k columns. We choose 
U to he a k X k unitary matrix. This choice has a considerable advantage. Let us recall 
the character formula [38] 

det( Xh^~'^~^^°') 

TtrX = —— — - — ^j^— — for X = diag(xi, ...,Xfc) and arbitrary R, (3.13) 

where Ra is the number of boxes in the a-th row of R. It implies that 

TtrtU = (detU)'' (3.14) 
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for the rectangular diagram R and for the k x k matrix By substituting fl3.12p and 
fl3.14p into (12.51) ■ we conclude that 

{WR)uiN) = ^ y^Mrff/ exp (^-^TrM^^ det(l + e*^^®[/"^)(dett/)". (3.15) 

By diagonalizing the matrices as 

M = diag(m,)^i, U = diag(e"iti, (3-16) 
and redefining U —U, the vev can be written as 

2N ^ 



{Wr)u{n) oc / JJ c/m„ JJ rfmi exp 

a i 

+ 5^1og (2sinh^^y + 5:io 



— ^ + ^ log(mi - TJij^ + n ^ Ma 

i=l i<j a=l 



^rrii—Ua 



(3.17) 



The constant of proportionality is independent of k and n, is negligible in the precision 
we work with, and is thus dropped. Since the integrand is analytic, the contours of 
integration may be deformed. Each Ua is integrated over a period 27ri. 

The leading behavior of the Wilson loop vev in the large limit can be computed 
in the saddle point approximation. Note that all the terms in the exponent in (13.171) 
are of order A^^. We thus expect a back-reaction of the m-eigenvalue distributions to 
the M-eigenvalues. This is in contrast to the case with a single row or column [13]. The 
saddle point equations are 

-^-^ + 2E^ E^7^^ = 0' (3-18) 

n^y coth ^^^^^ + y i = 0. (3.19) 

b^a i 

One can get useful intuitions by interpreting these equations as force balance conditions. 
The first term in (13.181) is a restorative force on mi. The first term in (I3.19P represents 
a constant force applied to each Ua- These are the external forces acting on the system 
of eigenvalues. The other terms are mutual repulsive forces among mj's and Ma's. 

In the absence of w-eigenvalues, the m-eigenvalues obey Wigner's semicircle law and 
spread over an interval of width 2\/\. This motivates us to assume that the m- and 



^ This can also be understood as follows. A column of length k gives the rank-fc anti-symmetric 
representation A^. The rectangular diagram R is obtained by symmetrizing n copies of A^. The 
representation Ak for U{k) is one dimensional and the trace in Ak is the determinant. Thus Tr^jr?/ = 
Trsym"(A.)t^ = T^A.»"C/ = (det C/)". 
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M-eigenvalues spread over regions of length scale a/A. Under this assumption that will 
be justified a posteriori, for large values of A, we can approximate various expressions as 
follows: 

-1 iiua<mi (3 20) 

if mj < Ma ' y ■ ) 



coth ~ = i -] if""<"V (3.21) 

2 1 +1 if Mfe < Ma 

We now propose an eigenvalue distribution that solves the saddle point equations. 
Suppose that the eigenvalues rrii split into the first group {m-^^|z = 1, ...,n} on the right 
and the second one {mp^|i = 1, N — n} on the left. All the Ma's are between them. We 
also assume that the two groups are far enough from each other so that l/{mf^^ — m^p) 
can be ignored. After we apply (13.201) and (13.211) . (I3.18P becomes 

- — mf ) + 2 J2 -m m + ^ = 0' (3-22) 

^^mF + 2 (,) ^ = 0. (3.23) 



A ""i I ^ (2) (2) 

Then the two groups individually obey Wigner's semi-circle law. The first group is spread 
over an interval of width l^fg^ = 0{VX) centered at kX/AN. For the second group the 
interval is centered at the origin with width 2^gs{N ~n) = 0{\^). fl319l) simplifies to 

^coth^^^^^ = 0. (3.24) 

This equation is solved by UaS uniformly distributed along a line in the imaginary direc- 
tion. The precise location of the line cannot be determined in this approximation. See 
Figure HI 

The parallel matrix model problem in the Chern-Simons case admits an exact solution 
for the finite 't Hooft parameter [30]. It exhibits an essentially identical distribution of 
eigenvalues. This gives us confidence in our proposal of the distribution. 

We now use the eigenvalue distribution to evaluate the Wilson loop vev. The exponent 
in (I3.17P can be evaluated, with the help of fl3.20p and fl3.2ip . as 

-X + ^ E + E - -5')' - X E(-f 

i i i<j i 

+ log(mf ) - mff + 0{N^ log A) 

i<j 

= ^ + 0{NHogX). (3.25) 
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N -n 



n 



Figure 4: The eigenvalue distributions in the rectangular case. The m-eigenvalues shown 
as black lines split into two groups, {rn^l^} on the right and {rrLp} on the left. The u- 
eigenvalues are distributed uniformly along the red line in the imaginary direction. 



The last term arises from completion of squares: 

- {2N/\){mf^Y + kmf^ = -(2iV/A)(mf^ - kX/4Nf + k^X/8N. 
We thus conclude that 

{WR)uiN) = exp i^-^ + 0{NHogX)j . (3.26) 

In the SU{N) case, we need to work with M' defined in (12.81) . Because of (12.91) . the 
Wilson loop vev is given by 

{Wr)su{n) = ^ J dMdUe-'-^^^' (det e^^)"^ det(l + ® f/-^)(det f/)". (3.27) 



(I3.25P is replaced by 



i<j 



The leading parts in A comes from completion of squares: 



A ' ' N"~' A V'^' 4:N'~ V NJ J ' 8N \^ N 



A'*' iV* AV 4AriV 

We find that 



n ^ 




AP 






SiV 


2 


AP 




+ 


8iV 


(^) 



(»H>..,., = -p("§(i-^)^('v-")i^(^)^o('v^'°sA)) 

= exp(^n(iV-n) + 0(iVMogA)J. (3.29) 



This is the special case of (11 ■2p . and is invariant under i? ^ i? or n ^ — n as it should 
be. 

Recall that each column in R represents a D5-brane. Since we have as many u- 
eigenvalues as columns, we identify the eigenvalues with D5-branes. The back-reaction 
of the m-eigenvalues to the existence of w-eigenvalues then represents the back-reaction 
of the geometry to the branes, i.e., geometric transition. Without operator insertions, 
there is a single black region representing AdS^ x S^. The hole in the region, which was 
identified with D5-branes in [T3], is created by the repulsion of and Ua in the present 
formulation. This is the Wilson loop analog of giant graviton branes represented by a hole 
in the fermion droplet [7j. In the topological and non-critical string literature [391 HO]. 
it is well-known that the insertion of a determinant operator similar to det(l -|- ® U) 
represents non-compact D-branes. In those cases the eigenvalues of U are the moduli of 
branes. Here there is a difference: the eigenvalues are first integrated over and are fixed 
only by saddle point equations. 



3.2 Eigenvalue bound states as D3-branes 

If we make use of the identity [371 EH] 

Tr/jV^Tr/je 



1 



R 



detll - 1/ (g) e 



instead of (I3.10p . we have another expression for Tr/je*^: 

1 



Tr^e 



M 



dV- 



-TtrV. 



(3.30) 



(3.31) 



det(l -y-i ®e^0 

dV is the Haar measure on the unitary group. More precisely, the contours need to be 
deformed so that the eigenvalues of V are larger than those of e^^ in magnitude. We 
choose to be an n X n unitary matrix, so that 

TtrV = (detV)'' (3.32) 

for the rectangular Young diagram R. We find that the Wilson loop vev is given by 
another 2-matrix model 



1 



dMdV exp 



{Wr)u{N) 

By diagonalizing the matrices, we get 
{Wr)u{n) oc j YldvaYldmiexp 



2N 



TrM^ 



1 



det(l - V-i (g)e 



M' 



-(detl/)^ (3.33) 



2N 

X 



N 



^ + ^ log(mi - mjY + k^Va 



1=1 



i<j 



a=l 



+ 5^ log ('2sinh^^y-^log(l 

a<b ^ ^ a,i 



(3.34) 
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Let us now analyze this model in the limit N ^ oo with X,n/N, and k/N finite. 
The saddle point equations are 

-— mi + 2> +> - = 0, 3.35 



and 



'3 



k + y coth "-^^^ - y = 0. (3.36) 

b^a i 

Again these equations can be interpreted as force balance conditions for the eigenvalues. 
The difference from the previous subsection is that the interaction between rrti and Va is 
attractive. 

Let us study the eigenvalue distribution for finite but large A. The ^-eigenvalues as 
a whole are pulled to the right by a force of magnitude kn. Unlike in the previous 
subsection, we cannot apply (13.201) and (I3.2ip to the terms involving mj and Va- If we 
did, Va would only be pulled to the right by m^, and the force balance would not be 
achieved. This suggests that each Va must be very close to mj. 

Each w-eigenvalue is pulled to the right by a force of magnitude k. To balance this, we 
assume that the f-eigenvalue and an m-eigenvalue sit very close to each other. We also 
assume, as will be justified a posteriori, that this distance is much smaller than all other 
length scales in the problem. There are n ^-eigenvalues and they are paired up with as 
many m-eigenvalues mi^ {i = l,...,n), forming n hound states. The bound states are 
denoted by {m-v)a, a = 1, ...,n. If we order the eigenvalues so that m^i^ ^ f i < ^2^^ ^ 
V2 < ... < TUn ^ Vn, fa IS pushcd to the right by vi, fa-i0 It is also pushed to the left 
by fa+i, I'n, but these forces are cancelled by the pull of rrin^ . Thus (13.361) in 

the approximation (13.211) becomes 

k + a ^—j-r + 0{N/\^) = 0. (3.37) 

Va-my 

The size of the a-th bound state is thus l/{k + a). 

Assuming that the group of bound states and the group of remaining m-eigenvalues 
{m- \i = 1, N—n} are far enough from each other, we can replace (13.351) by (I3.22p and 
(I3.23p . Since the m-eigenvalues are governed by the same equations as in the previous 
subsection, their distribution is the same. The distribution of the u-eigenvalues is identical 
to the distribution of m\^^ on the macroscopic scale. This then justifies the assumptions 
made in the argument. The computation of the Wilson loop vev is also identical. 



" Strictly speaking, for b such that < a-b < N/y/X ([3?2T1) cannot be apphcd. The error from ignoring 
the effects, however, is canceUed by the error from ignoring the forces from vt with < 6 — a < N/ \/X. 
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~ l/{k + a) 
•* ► 



Figure 5: The eigenvalue bound states {m-v)a- The distance between the two eigenvalues 
in the a-th bound state from the left is l/{k + a). The distance between neighboring 
bound states is much larger and is of the order \fX/N . 



n 

'- {{m-v)a} 

Figure 6: The eigenvalue distributions for M = diag(mj)^^ and V = diag(e''")^^^. {rrii} 
split into two groups {m[^^} and {ml^^}. The n w^'s are paired with as many m[^^'s to 
form n bound states {m-v)a- These bound states are distributed according to Wigner's 
semicircle law, as represented by the coincident black and red lines on the right. The 

(2) 

distribution of the remaining N — n m\ s independently obeys the semicircle law and is 
shown on the left. 



N -n 



Each row in the Young diagram represents a D3-brane [HI [15]. The rectangular 
diagram R then corresponds to n D3-branes. Since we have the same number of bound 
states, it is natural to identify each bound state with a D3-brane. This is the analog of 
an extra droplet in [7] where it was interpreted as dual giant graviton branes wrapping 
in AdS^ . 

4 General Young diagram 

This section deals with a general Young diagram of the form shown in Figure [2], with 
edge lengths rtj and kj all of order A^. We will obtain two multi-matrix models whose 
partition functions are the Wilson loop vev. 

4.1 Eigenvalue bound states as D5-branes 

Let us first see how to generalize the trick that led to 2-matrix models in the rectangle 
case. We begin with the expression 

Tr^e^^ = f dU^^^ det(l + e*'" ® U^^^~^)Tii^tUi. (4.38) 



12 



"J 







R = 



□ 

R^3) 



Figure 7: A shrinking sequence of Young diagrams R = R^^^ D R^'^^ D ... D R^^\ 



Here U^^^ G U{Ki). Recall that Kj = Yfj=ikj- dSH in this case allows us to write 

Tr^Tf/W = (detf/W)"iTr^(2)Tt/(^\ (4.39) 

Here -R*-^-* is obtained by removing the first ni rows from R. See Figure [71 Similarly, 

1 



Tr 



.f/(l) 



det(l -f/« ® t/(2)-i) 



(4.40) 



with f/(2) in [/(/sTs). To be more precise, to avoid singularities the integral is performed 
along the contours such that the eigenvalues of t/*^^-* are larger than those of U^^^ in 
magnitude. This time we have the relation 



Tr^(2)Tf/(2) = (detf/(2))«2Tr^(3)^f/(2). 



(4.41) 



By removing the first n2 rows from i?*^^-' one obtains R^^\ We now repeat the procedure 
as many times as we can. This yields 



Tr«e^ 

rff/(^' det(l + e^'^ ® f/(^)-^)(det U'^^^Y' I dU^^^ 



det(l -f/« ®m2)-i) 

dU^'^ -T^. TtAv. TTT^ (det f/^^^ )"« . 

det(l - f/(9-i) (g) f/(9)-i)^ ' 



^(detf/(2); 



"2 



(4.42) 



Here f/^^^ G U{Kj), and the integration contours are deformed so that the eigenvalues of 
U^^^ have larger absolute values than those of U^^~^\ 
We thus have 



X 



det(l + e*-'®f/«-^) 



det(l - t/W O f/(2)"i) ■■■det(l - f/(9-i) ® t/(9)-i) 



-. (4.43) 
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After redefining U^^^ —U^^^ and diagonalizing tlie matrices as M = diagi 
U^^^ = diag(e"'' )a=i^ ^^is becomes 



m. 



{Wr) oc J Yldm^lldui'^ exp 



2N 



N 



g Ki 



(I) 



i=l 



1=1 a=l 



log(mj — rrijY + log ( 2 sinli 

i<j I a<b 



log(i _ » ) - ^ 



/ a,b 



(4.44) 



Our aim is to understand the behavior in the hmit N ^ oo with X,nj/N, kj/N finite, 
for large values of A. The saddle point equations following from the action are 



AN 



A 



rii 



rii 



J^coth^^^^-^ 



1 

a=l '-^ 

^ 1 

■^^ — ' -ul, -rrii 



rrii — m,- ' — - ' -m.: 

a=l 



for z = l,...,Ar, 



i=l 



(2) (1) 



1 ■^T eV'-«^"' - 1 



for a 



coth — 



(/) (/) ^/-i 



y L_ 

/ V (/) (/-I) 

f ^ .111 ' —ti: ' 



1 



1 



(-f+l) (-T) 



(4.45) 



(4.46) 







for / = 2, ^ - 1 and a = 1, i^'/, (4.47) 



and 



n, 



9 + XI coth — 



(9) .,(9) 



for a = (4.4^ 



All the terms can be interpreted as forces on eigenvalues. Note that the interaction 
between and u"a^ is repulsive, while the one between xlI^ and is attractive. This 

suggests that some eigenvalues form bound states. 

Let us work in the A ^ 00 approximation. The following distribution of eigenvalues 
solves the saddle point equations. It involves g types of bound statesjfl The J-th type 
of bound state, which we denote by {u^^^-u^'^^-...-v!^^^\ contains a single u^'^^-eigenvalue 
for 1 < J < /, as shown in Figure [81 The m-eigenvalues split into + 1 groups 
{rti\^\i = The kj (M'^-^)-M*^^)-...-M('^))-bound states are distributed along a line 



Though it is convenient to talk about bound states, the size of a bound state is in fact of the same 
order 0{1/N) as the distance to the neighboring bound state. 



14 




Figure 8: There are g types of bound states {u^^^-u^'^^-...-u'^^^), 1=1, ■..,g- 



in the imaginary direction, and separate {Tn[^^} from {m[^~^^^}. See Figure [91 For inter- 
actions between two bound states or between a bound state and an m-eigenvalue, we can 
apply (13.201) and fl3.2ip . The (u(^^-u^^^-...-n(^^)-bound state is pulled to the right by an 
external force of magnitude ni+n2 + --.+ni, and pushed back to the left by ni+n2 + --.+ni 
m-eigenvalues from {ml'^^} for 1 < J < I . Bound states of different kinds do not in- 
teract with each other: forces cancel out among constituent eigenvalues. The saddle 
point equations are satisfied if the bound states are uniformly distributed along verti- 
cal lines. Calculation of relative positions of constituent eigenvalues in a bound state is 
straightforward. In the bound state {u^^^-u^'^^-...-u^^^), the distance between the u^"^^- and 
„(^+i)_eigenvalues is 1/(EJ^=j+i + Ejfij ^i^) = 0{1/N) for 1 < J < / - 1. 



Ug+l 



n„ 



n2 



ni 



Figure 9: The eigenvalue distributions for M, [/(i), [/(^-i), and U'-^l The N 



m-eigenvalues split into 5^ + 1 groups {m. 



1, n/}, I = 1, g + 1. The kj {u 



u 



(2). 



-U 



-bound states are distributed uniformly between {mf^} and {mf'^^^}, and 



are represented by / coincident colored lines in the imaginary direction. 



It seems reasonable to believe, based on the analogy with [30], that the proposed 
eigenvalue distribution is the leading saddle point. Let us calculate the Wilson loop vev 
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by using the proposed distribution. The exponent of (14.441) becomes 

^ g+l nj g nj g+1 

a" 



1=1 i=l 1=1 1=1 1=1 i<j 

9 



J2^njK] + 0{NHogX). (4.49) 



8N 
1=1 



We get the final result for the U (N) Yang-Mills: 

In the SU{N) case, (12. 9p allows us to write 

{V^R)suiN) = ^y"c/Mrff/(i)rft/(2)...rff/(9)e-^TrM^(dete^ 



,M^-\B\/N 



(4.51) 



X (det [/(^))"ndet f/(2))"2...(det f/^^^)"'' det(l + ® f/^^^"^ 

1 1 
^ det(l - f/« ® f/(2)-i)'"det(l - ® f/(9)-i) ■ 

(I4.49P is replaced by 

/=1 i 1=1 i 

9+1 

+ E E log(mf ) - ))2 + 0(iV2 log A) 

7=1 i<j 
9+1 . 

= E^^^(^^-|^l/^)' + ^(^'l°g^)- (4.52) 

7=1 

Recall that Kg^i = 0. This is the result (11.21) presented in the introduction. 

4.2 Eigenvalue bound states as D3-branes 

Another matrix model can be obtained by starting with the formula 

Tr«e*^= fdV^'^—- ^L^^^Tr^^W. (4.53) 

J det(l - e*^ ® yW-i) ^ ^ 

Here V^^^ G U{Ni) and dV^^'' is the Haar measure. We have defined Nj := Ylj~=\~^ ^J- 
Note that 

TtrV^^') = (detl^(^))^'''TrQ(2)y(^\ (4.54) 
16 



iVi 



No 



R = Q('^ g(2) . . . Qig-l) Qig) 

Figure 10: Another shrinking sequence of Young diagrams R = Q^^^ D Q^"^^ D ... D Q^^^. 



Here Q^'^^ is obtained by removing the first kg rows from R. See Figure [TOl Similarly, 



det(l-V^W 1^(2) 



(4.55) 



with V^"^^ essentially in U{N2). To avoid singularities the integral is performed along 
appropriate contours. We repeat the same procedure as many times as possible. This 
gives 



Tr/iC 



M 



det(l-e^^®V«-i)^ ^ J det(l - yd) ® V^(2)-i) 



X 



(detV^(2))fc 



s-i 



det(l - V^9-i) (g) 



-(detl^(^'))'=i. (4.56) 



Here V^'^'^) G U{Nj), and the integration contours are deformed so that the eigenvalues of 
V^^^ have larger absolute values than V^^~^K 
We thus have 

{Wr) = ^ J dMdV^^UV^^\..dV^'^e-^^''''\det l^«)'=«(det V^^'>)''^-\..{det V'^^^)''' 

^ det(l - ® F(i)"i) det(l - ® - det(l - V(a-^) ® V^(^^)-i) ' ^^'^'^^ 

In terms of the eigenvalues, the Wilson loop is given by 



2N 



N 



g Ni 



in 



i=l 



1=1 a=l 



+ ^2 ^°s("^i - 



(n jn' 



Va -V 



i<j 



I a<b 



i,a / a, 6 



(4.58) 
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Let us study this model in the large limit with nj/N and kj/N fixed. The saddle 
point equations are 



and 



AN 2 1 

+ y, + y, — = for ^ = 1, N, (4.59) 

b^a i=l ^ fe=l e (J i 

for a = l,...,7Vi, (4.60) 

K+i-i + ^coth " ^ ^ - 5Z ^ + E ,(^+1) ,(^) , = 

for / = 2, ^ - 1 and a = 1, Nj, (4.61) 



(a) _ (9) ^s-i 1 
/ci + 5^ coth - 5^ = for a = 1, iV„ (4.62) 

b^a ^ b=l -1 

The forces between rrii and v^a'^ as well as between Va^ and w^''^^"^^ are attractive. 




Figure 11: There are g + I types of bound states {m-v^^^-...-v^^^), / = 0, 1, ...,g. 



Our proposed distribution of the eigenvalues for large A is as follows. There are g types 
of bound states denoted by {m-v^^^-...-v^^^) for J = 0, 1, ■■.,g. Each bound state of the 
/-th type has a single m-eigenvalue as well as a single u*^ ■^-'-eigenvalue for 1 < J < I . The 
size of a bound state is of the order 1/A^, and is much smaller than the distance between 
two neighboring bound states. The Ug+i-i (m-t;*^^)-...-w^'^^)-bound states are distributed 
according to the semi-circle law centered at XKj/AN. Using fl3.20p and (13.211) . one can 
confirm that the force balance among bound states is achieved. It is straightforward 
to calculate the distance between eigenvalues in a bound state. In the a-th (m-v^^^- 
...-v^^^) bound state, the distance between the v^"^^- and f '^'^"'"^^-eigenvalues turns out to 
be l/((/ - J)a + J2K=jh+i~K) = 0{1/N) for < J < / - 1. Here = m and 
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kg+i = 0. This configuration solves the saddle point equations. The distribution of the 
m-eigenvalues is identical to the one in the previous subsection, and leads to the same 
results for the Wilson loop vev. 



Figure 12: The eigenvalue distributions for M, V^'^\ V^^\..., V^3-^\ V^^\ The Ug+i-i 
(m-f *^^)-...-f *^'^^)-bound states are distributed according to Wigner's semi-circle law. 
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